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In f 2 -based heavy fermion systems with a crystalline-electric- field (CEF) singlet ground state, 
the Non-Fermi Liquid (NFL) arises around the quantum critical point (QCP) due to the com- 
petition between the CEF singlet and the Kondo-Yosida singlet states. In such a case, the 
characteristic temperature Tp at which the entropy starts to decrease toward zero is sup- 
pressed by the effect of the competition, compared to both energy scales characterizing each 
singlet state, the lower Kondo temperature (Tk2) and the CEF splitting (A). We show that 
in the case of tetragonal symmetry Tp is not affected by the magnetic field up to Hi which is 
determined by the distance from the QCP or characteristic energy scales of each singlet states, 
not by Tp itself. As a result, in the vicinity of QCP, there are parameter regions where the 
NFL is robust against the magnetic field, at an observable temperature range T > Tp, up to 
if* which is far larger than Tp and less than min(Tk2, A). Our result suggests that such an 
anomalous NFL behavior can arise also in systems with other CEF symmetry, which might 
provide us with the basis to understand the anomalous behaviors of UBei3. 

KEYWORDS: numerical renormalization group, heavy fermion, crystalline-electric-field singlet, f 2 , non- 
Fermi liquid, Kondo effect, magnetic field effect 



1. Introduction 

In the last decade or so, non-Fermi liquid (NFL) be- 
haviors around quantum critical point (QCP) have been 
one of main issues in physics, not only in heavy fermion 
systems, 1 ) but also in those exhibiting the Mott transi- 
tion. 2 ) Of these NFL behaviors, those of heavy fermion 
systems with ^-configuration form a kind of subclass in 
which the QCP is triggered by local criticalities: such 
as the two-channel Kondo effect (TCKE) due to the 
non-Kramers doublet state, 3 ' 4 ) and that caused by the 
competition between the crystalline-electric field (CEF) 
singlet and the Kondo-Yosida (K-Y) singlet states. 5, 6 ) 
The former TCKE was reported to be observed in 
Lai-^Pr^Pba that has a T3 non-Kramers doublet ground 
state in the cubic symmetry. 7 ) The NFL behaviors in 
Thi_ :E U x Ru2Si2 were understood in a unified way by as- 
suming that the system is located near the phase bound- 
ary between the CEF singlet and the K-Y singlet states. 5 ) 
owever, a detailed study about the magnetic field depen- 
dence on NFL behaviors has not been performed so far. 

In the present paper, we investigate the magnetic 
field dependence of NFL behaviors in the specific heat 
Cimp(^) and the entropy S'imp(^) due to f-electrons with 
the two-orbital impurity Anderson model in a tetrag- 
onal symmetry with the CEF singlet ground state on 
the basis of the numerical renormalization group (NRG) 
method. 8,9 ^ We discuss how the magnetic field, H z , 
changes the characteristic temperature, Tp, which is 
defined as the temperature at which the temperature 
derivative of entropy, 95 , i mp (T)/9(log T), takes the max- 
imum value as 5'i m p(^) approaching as T — >• 0. In the 



vicinity of the QCP, Tp is suppressed by the effect of 
the competition between the CEF singlet and the K-Y 
singlet states for H z = 0, and the NFL behaviors occur 
at Tp < T < Tk2, where Tk2 is the lower Kondo tem- 
perature of two orbitals, as in the case of TCKE. The 
magnetic field is shown not to affect Tp up to a certain 
value HI which is determined approximately by the con- 
dition that the effect of the magnetic field, destroying a 
criticality of the TCKE type, becomes comparable to the 
effect of the deviation from the criticality at H z = 0. JT* 
so determined is far larger than T^(H Z = 0) for a reason- 
able set of parameters. As a result, the NFL behaviors 
become robust against the magnetic field up to if* ~ Tk2 
which is about hundred times larger than T^(H Z = 0). 

This paper is organized as follows. In §2, the model 
Hamiltonian is introduced and transformed into a form 
suitable for the NRG calculation. In §3, we discuss how 
the characteristic temperature Tp is affected by the ef- 
fect of the competition between the CEF singlet and the 
K-Y singlet states in the case of H z = 0. In §4, we 
demonstrate the magnetic field dependence of Tp and 
7im P (^) = Ci mp /T. In the vicinity of the QCP, there 
are parameter regions where — log T behavior of 7i mp , at 
temperature Tp < T < min(Tx, A), is robust against the 
magnetic field. In §5, we investigate how such an anoma- 
lous NFL is affected by the change of the characteristic 
energy scale of two singlet states. In §6, we summarize 
our results and discuss their applicability for understand- 
ing the magnetically robust NFL behaviors observed in 
UBei3 because such an NFL being robust against the 
magnetic field can arise in systems with other symme- 
try if the K-Y singlet state and the CEF singlet state 
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compete for the ground state. 

2. Model Hamiltonian 

In this section, we recapitulate discussions of ref.5 
about how to derive the model Hamiltonian for dis- 
cussing the competition between the K-Y singlet and the 
CEF singlet states in f 2 -configuration on the basis of the 
j — j coupling scheme in the tetragonal symmetry. We 
restrict the f 1 state within two low- lying doublet states 
out of three doublets of j = 5/2 orbitals, and allot the 
pseudospin representation for these states as follows: 
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Here, for example, ^,0) represents the state where or- 

(2) 

bital 1 (Ty ) with up pseudospin is occupied and orbital 
2 (Tq) is empty. We also restrict the f 2 state within four 
low-lying states out of states allowed in J = 4 manifold, 
and construct these four states with the direct product 
of f 1 states. Here, we have discarded states where two 
f-electrons occupy the same orbital, |ti?0)? |0>ti)> be- 
cause the intra-orbital Coulomb repulsion is larger than 
the inter-orbital one. Then, low-lying four f 2 states are 
expressed as 

|r 4 > = 



|r 3 > 

Hfl> 
|if>> 



7T 



2)-|-2» = ^(|+,t)-|t,l», (5) 

2> + |-2» = ^(|t,|) + |l,t», (6) 

/3|+3>-a|-l) = |t,t>, (7) 

f3\-3)-a\ + l) = \i,i). (8) 



It is noted that we cannot determined coefficients, a and 
/?, because we have discarded one of the doublet in f 1 - 
configuration. Therefore, in this paper, we take its j — j 
coupling representation as f 2 states with symmetry 
as shown in Appendix including the derivation of eqs.(l)- 
(8). 

We assume that the CEF ground state is the singlet 
(T4), the first excited CEF states are magnetic doublet 
(r 5 ) with the excitation energy A, and the second ex- 
cited CEF state is the singlet ^3) with the excitation 
energy as shown in Fig.l. Such a CEF level scheme 
can be reproduced by introducing the "antiferromagnetic 
Hund's-rule coupling" for the pseudospin as 

2 



H Hu nd = -± [SfSz + SfS+J + J Z SISI (9) 



where coupling constants are defined as J± = K and 
J z = 2A — K, respectively, and Si is a pseudospin op- 
erator of the localized electron in the orbital i defined 
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Fig. 1. CEF level scheme of low- lying f states and their eigen- 

states. 
as 
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Furthermore, assuming that f-electrons constructing the 
f 2 state hybridize with conduction electrons which have 
the same symmetry as each f 1 state. Thus the system can 
be described by the two-orbital impurity Anderson model 
with the "antiferromagnetic Hund's-rule coupling" as fol- 
lows: 



Ti — He + Hhyb + Tif + HHund, 



(11) 

(12) 
(13) 



% = E Y. E ^l^° + E E y (14) 

i=l,2 cr i=l,2 cr 

where fi a (fi a ) and c^ ia (cl.^) are annihilation (creation) 
operators of the f-electron on the orbital i with the en- 
ergy Efi and the conduction electron with wave vector 
k hybridizing with the f-electron with the symmetry of 
the orbital i with strength V^. Here, the on-site intra- 
orbital Coulomb repulsion Ui is explicitly taken into 
account, while other Coulomb repulsion terms like the 
inter-orbital or the exchange interaction, are implicitly 
included in the "antiferromagnetic Hund's-rule coupling" 
of (9). 

To analyze properties of the system described by the 
Hamiltonian (11) by the Wilson NRG method, 8,9 ) we 
transform the conduction electron part as usual. For sim- 
plicity, we take conduction bands to be isotropic in mo- 
mentum space, i.e. the hybridization depends only on the 
orbital z, = Vi, and symmetric in the energy space 
(with an extent from — D to D) about the Fermi level. We 
discretize conduction bands logarithmically with the dis- 
cretization parameter, A, and perform the unitary trans- 
formation assuming the density of state in conduction 
bands as constant. Thus, eqs. (12) and (13) can be rewrit- 
ten as 

00 



(16) 
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where (fj n ) is the annihilation (creation) operator 
of the conduction electron in the shell orbital whose ex- 
tent is /cpA n / 2 and ji-\ a = fia- The hopping integral 
between n-th and (n + l)-th shell states, t n , is expressed 
as 

L>(l + A- 1 )(l-A- n - 1 ) 



t r . 



2y/(l - A- 2 ™- 1 )(l - A" 2 "" 3 ) ' 



(17) 



Then, we define Hn which approaches %/(D(l + A 1 )/2) 
in the limit N — >• oo as follows: 



: A (iV-l)/2 



IV — i 



. (18) 



where the tilde indicates that energies are measured in 
a unit of D(l + A _1 )/2. The Hamiltonian (18) satisfies 
the recursion relation 

(19) 

We solve the whole sequence of Hamiltonian (Hat) by 
using the recursive form (19) with keeping states up to 
1500 states in each iteration step, and use A = 3.0 in all 
the calculations below unless explicitly stated. 

3. NFL Behavior due to Competition between 
CEF and K-Y singlets 

In this section, we discuss the effect of the competi- 
tion between the CEF singlet and the K-Y singlet states, 
which can give rise to a NFL state. It is already known 
that the system described by the Hamiltonian (11) has 
the competition between the K-Y singlet and the f 2 - 
CEF singlet states. 5 ) In general, the energy level and 
the strength of hybridization with conduction electron 
in each f-orbital are different. In the present paper, we 
take parameters so that the Kondo temperature of or- 
bital 2 is always lower than that of orbital 1: i.e., we set 
parameters of the two-orbital impurity Anderson model, 
eq (11), as E fl = E f2 = -0.4, Ui = U 2 = 1.0, V x = 0.45, 
and V2 = 0.3. Hereafter, the unit of energy is taken as D 
unless stated explicitly. In the case of K = A = 0, the 
model Hamiltonian, eq. (11), reduces to two independent 
impurity Anderson models. The Kondo temperature of 
each orbital can be determined by the Wilson's defini- 
tion, 4TKXim P (^ = 0) = 0.413, for conventional Ander- 
son model as T K1 = 6.10 x 10~ 2 and T K2 = 6.01 x 10~ 3 , 
respectively. 

For the finite value of CEF parameters, (K, A), there 
are two stable Fermi Liquid (FL) fixed points correspond- 
ing to two singlet ground states as shown in Fig. 2: the 
K-Y singlet (filled circles) and the CEF singlet (open cir- 
cles) fixed points. At the boundary of these two regions 
of FL fixed points, there exists a curve of critical points, 
across which energy spectra for even and odd iteration 
interchange, and NFL behaviors appear in the vicinity 
of the boundary. To analyze further, we fix one of the 
CEF parameters as K = 0.16, and calculate the physi- 
cal properties for a series of values of the CEF splitting 
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Fig. 2. Phase diagram of the ground state in K — A plane. Filled 
circles represent the K-Y singlet fixed point and open circles 
represent the / 2 -CEF singlet fixed point. Parameter set is Ef 1 = 
E f2 = -0.4, U!=U 2 = 1.0, Vi = 0.45, and V 2 = 0.3. 
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Fig. 3. Temperature dependence of the entropy due to f-electrons 
in systems near the critical point. Parameter set is the same as 
that used in Fig. 2. In order to obtain the result with a higher 
accuracy, 3000 states are kept in each step of NRG. Ground 
states of each system are indicated by open symbols for the K-Y 
singlet, and filled symbols for the CEF singlet. The characteris- 
tic temperature Tp is given by that making dSi mp (T) / d(\og T) 
maximum at the lower temperature side. 

parameter A. Analyzing near the critical point in more 
detail, the critical value of A is determined as A* ~ 0.112 
for K = 0.16. 

Fig. 3 shows the result of the temperature dependence 
of Simp CO > the entropy due to f-electrons, near the criti- 
cal point. As mentioned above, the characteristic temper- 
ature Tp is defined as the temperature at which the tem- 
perature derivative of entropy, <9Si mp (T)/<9(log T), takes 
the maximum value just before Si mp (T) approaching 
as T — ^ 0. As seen in Fig. 3, Tp is drastically suppressed 
by the effect of the competition near the critical value of 
CEF splitting A = 0.112 ~ A*. 

Fig. 4 shows the A dependence of Tp which is ob- 
tained by numerical calculations of Si mp (T). In the case 
of A < A*, the K-Y singlet state is the ground state, and 
two localized moments Si and S2 are screened out inde- 
pendently by corresponding conduction electrons, where 
each Kondo temperature is affected by the interaction 
between f-electrons. In this case, the total phase shift 
of conduction electrons characterizing this fixed point is 
5 = 7r (Si = 7r/2, S2 = tt/2), and Tp is given by a value 
slightly lower than the Kondo temperature Tk2, if A is 



4 J. Phys. Soc. Jpn. 



Full Paper 



Author Name 



10" 



10" 



10"' 



10 



i— 

r D \ 
: □ 




□ i 

□ 

r □ 

□ 

r \ 

: □ 

■ \ 


/ 

/ 
d 
□ 

/ 

□ 


■ K-Y Singlet \ 

■ i i i i 


/ CEF Singlet 

i i i 



K=0.16 



0.06 0.08 0.10 0.12 0.14 
A 

Fig. 4. A dependence of Tp. The effect of the competition be- 
tween two singlet states suppresses Tp , and in particular Tp = 
at the critical point A = 0.112 ~ A*. 




Fig. 5. Schematic energy levels of two singlet ground states. The 
CEF singlet state is stabilized relative to the K-Y singlet state 
as A increases. 

much smaller than A*. On the other hand, in the case of 
A > A*, the CEF splitting (antiferromagnetic interac- 
tion between f-electrons in the model Hamiltonian, (11)) 
is so large compared to the energy gain related to the 
formation of K-Y singlet states that the CEF singlet be- 
comes the ground state. In this case, the remaining con- 
duction electrons are not scattered by f-electrons, and as 
a result the total phase shift is 5 = (5i = 0,#2 = 0). 
When A ^> A*, Tp becomes close to the excitation en- 
ergy K between two singlet states. 

Such an interchange of the ground state can be un- 
derstood by considering that the increase of A causes 
the stabilization of the level of the CEF singlet state as 
shown in Fig. 5. In the case of A ~ A*, Tp is deter- 
mined not by characteristic energies of the K-Y singlet 
and the CEF singlet states, but by the energy splitting 
between two singlet states, AE: i.e., Tp ~ AE. Partic- 
ularly, at the critical point, the degeneracy of the K-Y 
singlet and the CEF singlet states is not lifted even at 
T = 0, making Tp = and lim T ^ Simp = 0.5 log 2. In 
other words, at low enough temperatures, the localized 
moment Si of orbital 1 has already been screened out by 
conduction electrons in orbital 1 below Tki, while §2 of 
orbital 2 still has the degree of freedom as localized mo- 
ment. Therefore, the effective Hamiltonian of (11) near 
the fixed point behaves as the two-channel Kondo model 
(TCKM) 10 ' 11 ) because S 2 interacts with two "conduc- 
tion" electron channels, one is the conduction electrons 
on orbital 2 and the other is a complex of conduction 
electrons on orbital 1 and screened Si as discussed in 
ref. 10. 

Fig. 6 shows the A dependence of the Sommerfeld co- 




Fig. 6. A dependence of the Sommerfeld coefficient 7i m p(T) = 
Ci mp (T)/T due to the f-electrons for various temperature. The 
ground state switch at A = A* ~ 0.112 from the K-Y singlet 
ground state for A < A* to the CEF singlet for A > A*. 

efficient, 7i mp (^) = C{ mp (T)/T, due to f-electrons for 
various temperatures. For all A shown in Fig. 6, 7i mp (^) 
increases monotonically down to T = 7.0 x 10 _T as de- 
creasing T. At A = A* ~ 0.112, the increase of 7i mp (T) 
does not stop and exhibits divergence in the limit T — »• 
because the structure of the fixed point is the same as 
that of TCKM as discussed above. For A off the crit- 
ical value A*, the increase of 7i mp (^) stops at around 
the characteristic temperature Tp leading to the Fermi 
liquid behavior at T < Tp. 7i mp (^) takes a dip struc- 
ture around A ~ A* at higher temperature region. This 
is because Si mp (T) has only a weak T dependence in a 
wide temperature range ~ Tp < T < Tk2 or A around 
A ~ A* as can be seen in Fig. 3. 

It is remarked that the enhanced part of 7i mp (T) near 
A ~ A* in the low temperature limit from the back- 
ground part at |A — A*| ^> A* arises from the effect of 
the competition between the K-Y singlet and the CEF 
singlet states. The part of the background is essentially 
given by an inverse of Tk2 or A, and is overwhelmed by 
the enhanced part near A ~ A*. Note that the ordinate 
of Fig. 6 is represented in a logarithmic scale. 

Although we take A as a control parameter here, we 
can expect a similar behavior of 7i mp through other pa- 
rameters, such as the hybridizations Vi and V2, which can 
also control the competition between levels of two singlet 
states. It is also remarked that such an anomalous behav- 
ior of 7i mp can be realized in systems with other symme- 
try: e.g., in UBei3 with cubic symmetry. 12 14 ) In this 
material, 7 shows the similar behavior as shown in Fig. 6 
through the change of the lattice constant, ao, which is 
controlled by replacing the U atom partly with other 
nonmagnetic elements. It is remarkable that 7 takes a 
maximum value at ao = ag, which is approximately the 
same as the lattice constant of UBei3. 15 ) Experimentally, 
in a series of materials with ao < a^, the Kondo like up- 
turn is observed in the resistivity in the low temperature 
region, while in those with ao > a^, the temperature de- 
pendence of the resistivity can be explained by the effect 
of the CEF with the singlet ground state. Then, we ex- 
pect that UBei3 is located near the critical point in this 
series of materials. 
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Fig. 7. Magnetic field dependence of Tp near the critical point. 
The parameters related to f-electrons are the same as Fig. 2. Cir- 
cles indicate characteristic magnetic fields Hz*'s. 

4. Magnetic Field Dependence of Non-Fermi 
Liquid Behavior 

In this section, we discuss the magnetic field de- 
pendence of the NFL behavior of 7i m p (T). The effect 
of the magnetic field on f 1 states is taken into ac- 
count through the Zeeman term for total angular mo- 
ment, Hzeeman^ 1 ) = -gj^Bj z H z , with j = 5/2 and 
gj = 6/7. That on f 2 states arises from the diag- 

(2) 

onal (for T K 5 ' doublet) and the off-diagonal (for T3 
and T4 singlets) matrix elements of two-electron Zee- 
man term H Z eeman(/ 2 ); e.g., (T^ |%Zeeman(/ 2 ) \^±) = 

^llgjH B H z /7, and (r 3 |Hzeeman(/ 2 )|r 4 ) = -2gjH B H z . 

Here, %z eemant/ ) Consists of two ^Zeeman (/ )• 

In Fig. 7, we show the magnetic field dependence of the 
characteristic temperature Tp near the critical point; i.e., 
A = 0.108, 0.110, 0.112(~ A*), 0.114, 0.116 and 0.118. It 
is noted that T^(H Z ) remains constant for H z less than 
the characteristic magnetic field H* which is defined ap- 
proximately as that from which T^(H Z ) starts to increase 
as increasing H z (as shown by circles in Fig. 7). Explic- 
itly, the characteristic magnetic field if*'s are given as 
H* ~ 3 x 10~ 4 for A = 0.106 and 0.118, H * ~ 2 x 10~ 4 
for A = 0.108 and 0.116, H * ~ 3 x 10~ 5 for A = 0.110 
and 0.114, and H * ~ 1 x 10~ 5 for A = 0.112. H* has a 
tendency of approaching zero as the critical fixed point is 
approached, i.e., A — >> A*. For CEF parameter A shown 
in Fig. 7, H* is much smaller than the lower Kondo tem- 
perature Tk2 — 6.01 x 10 -3 , so that the magnetic field 
H z < H* has little influence on the K-Y singlet state. 
Then, H* is considered to be determined by a competi- 
tion of two effects which destroy the TCKM-type NFL 
fixed point: one is a distance of A from A* and the 
other is the magnetic field which breaks the degener- 
acy corresponding to *Si mp (T = 0) = 0.5 log 2 due to the 
TCKE, the origin of the TCKM-type NFL fixed point. 
Namely, H* is given by the energy scale characterizing 
a crossover from the TCKM-type NFL behavior to the 
polarized Fermi liquid behavior beyond the effect of the 
distance of A from the critical value A*. Since 7i m p(T) 
exhibits the divergent increase around A ~ A* in the 
temperature region T > T^(H Z ) as decreasing T, 7i mp (T) 
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Fig. 8. Temperature dependence of 7 for (a) A = 0.112 (~ A*) 
and (b) A = 0.118 under various magnetic fields. In the case of 
(b), the NFL behavior of 7 is robust against a magnetic field of 
up to H z = 1.2 x IO -3 for T > 3.0 x 10~ 5 in spite of Tp ~ 
1.69 x 10 -5 . The parameters related to f-electrons are the same 
as those used in Fig. 2. 

exhibits a NFL behavior in the same temperature region 
T > Tp(H z ). Since Tp(H z ) remains almost unchanged up 
to H z = H*, the NFL behaviors are expected to remain 
robustly even under the magnetic field H z > T^(H Z ) 
so long as H z < Tk2- This behavior is reproduced by 
explicit calculations of 7i mp (^) under various magnetic 
fields as shown below. 

In Fig. 8, we show the temperature dependence of 
7im P (T) for A = 0.112 (~ A*) and A = 0.118 under var- 
ious magnetic fields of up to H z = 1.2 x 10 -3 . Extremely 
close to the criticality at A = 0.112 A*, 7i m p(^) is 
enhanced by the magnetic field as shown in Fig. 8 (a). 
This is because T^(H Z ) increases appreciably from 10 -7 
to 10 -5 corresponding to the increase of the magnetic 
field H z from 10 -4 to 10 -3 , resulting in an increase of 
0S imp (T)/0(log T) = Cimp(T), so 7 imp(T), at T > IO" 5 . 
On the other hand, at A = 0.118 slightly off the criti- 
cality, 7im P (^) is robust against the magnetic field up to 
H z = 1.2 x 10 -3 for the temperature region T > 3 x 10 -5 
as shown in Fig. 8(b). This is because T^(H Z ) remains al- 
most unchanged up to H z = H* ~ 10 -3 so that 7i mp (^) 
remains the same as that at H z = for T > Tp ~ 10 -5 . 

These kinds of NFL behaviors arise also in the region 
of the K-Y singlet state, i.e., A < A*, although we do 
not show the results explicitly. 

5. Kondo- Temperature Dependence of Non- 
Fermi Liquid Behavior under Magnetic Field 

In this section, we investigate the properties of the 
NFL behavior of 7i mp (T) under magnetic field of sys- 
tems with other Tk2 by changing V2 as V2 = 0.25 and 
0.20 for various sets of the CEF parameter, (if, A). Other 
parameters are set to be the same as those in the pre- 
vious section: i.e., Efi = Ef2 = — 0.4, JJ\ = U2 = 1.0, 
and V\ = 0.45. In the case of K = A = 0, each lower 
Kondo temperature can also be determined by the Wil- 
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Fig. 9. Temperature dependence of 7i m p(r) for a series of mag- 
netic fields in the system for the hybridization (a)V2 = 0.25 and 
(h)V 2 = 0.20. 

son's definition as T K2 = 1.27 x 1(T 3 for V 2 = 0.25 
and T K2 = 8.92 x 1(T 5 for V 2 = 0.20, respectively. To 
analyze further, we also fix one of the CEF parame- 
ters as K = 0.16 and calculate 7i mp (T) for a series 
of A under various magnetic fields. It is natural that 
A* (corresponding to the critical point) becomes small 
with decreasing Tk2 because the energy gain due to 
the formation of the K-Y singlet state decreases with 
a smaller V 2 . The critical value of A is determined as 
A* ~ 0.054 for V 2 = 0.25 and A* ~ 0.024 for V 2 = 0.20, 
respectively. Fig. 9 shows the temperature dependence 
of 7i mp (T) of the system with the CEF ground state: 

(a) A = 0.062 > A* ~ 0.054 for V 2 = 0.25 and (b) 
A = 0.032 > A* ~ 0.024 for V 2 = 0.20. The NFL be- 
havior being robust against the magnetic field occurs in 
a temperature region of T > Tp up to H z ~ H* in the 
former case (a), while in the latter case (b) the magnetic 
field has considerable influence on the NFL behavior. 

For H z = 0, Tp is also suppressed as in the case of 
V 2 = 0.30 in the vicinity of the critical point A ~ A*. It 
is noted that the decrease of Tk2 and A does not appre- 
ciably affect Tp, i.e. Tp ~ 10 -5 for both cases of (a) and 

(b) , which is determined from calculations correspond- 
ing to Fig. 3. This is because Tp is determined by the 
energy splitting between the K-Y singlet and the CEF 
singlet states, and does not depend on the characteristic 
energy scale of each singlet state. Under the magnetic 
field, the effect on the NFL behavior is markedly dif- 
ferent in two cases (a) and (b). In the case of (a) with 
V 2 = 0.25, the NFL behavior of 7i m p(T) is rather robust 
against the magnetic field (up to H*) in a wide temper- 
ature range (T > Tp) as in the case of V 2 = 0.30, while 
in the case of (b) with V 2 = 0.20, 7i mp (T) is sensitive to 
the magnetic field because the characteristic magnetic 
field H* is comparable to the lower Kondo temperature, 
Tk2- Namely, in the case of V 2 = 0.20, the magnetic field 
H z > T K 2 — 8-92 x 10 -5 suppresses 7i mp (T) by breaking 



the K-Y singlet ground state. It is noted that in the case 
Tk2 > A, the suppression of 7i m p(T) as in the case of 
Fig. 8(b) is expected for H z > A by breaking the CEF 
singlet states. Thus, the magnetic field dependence of 
the NFL behavior of 7i mp (T) is determined not by the 
characteristic temperature Tp, but by the characteristic 
magnetic field H* which is determined by the character- 
istic energy scale of each singlet state, Tk2 and A, or the 
distance from the critical point. 

6. Conclusion and Discussion 

We have investigated the effect of the magnetic field on 
the NFL behaviors due to the competition between the 
K-Y singlet and the CEF singlet states in f^-based heavy 
fermion systems with tetragonal symmetry. The effect of 
the competition suppresses the characteristic tempera- 
ture Tp, corresponding to a peak of the specific heat, 
Ci m p(r), to a much smaller value than the characteris- 
tic energy scale of each singlet states: i.e., Tk2, the lower 
Kondo temperature, and A, the energy splitting between 
the CEF singlet ground state and the first excited dou- 
blet states. Tp is determined approximately by AE, the 
energy difference between two singlet states, and there 
exists the two-channel Kondo model (TCKM) type NFL 
behaviors at Tp < T < Tk2- Namely, near the critical 
point, A ~ A*, the Sommerfeld coefficient 7i mp (T) ex- 
hibits a NFL behavior (7 imp (T) ex - log T) at T > T*. 

In the vicinity of the critical point, Tp was shown not 
to be affected by the magnetic field up to a certain value 
iJ*, while Tp is increased for H* < H z < min(T K 2, A). 
As a result, the NFL behavior of 7i mp at T > Tp is robust 
against the magnetic field H < H*. Then, for reason- 
able sets of parameters, the NFL behaviors being robust 
against a magnetic field of up to H* can occur at an 
observable temperature range. Thus, the magnetic field 
dependence of this NFL is characterized by H* which 
is determined by the characteristic energy scales of two 
singlet states and the distance from the critical point. 

In the present paper, we have discussed physical prop- 
erties in the tetragonal symmetry. However, also in the 
case of other crystal symmetries, it is expected that there 
remains the effect of the competition between the K-Y 
singlet and the CEF singlet states, leading to the NFL 
behaviors similar to the present case. One example would 
be the case of the cubic system UBei3 which seems to be 
located near the phase boundary between the K-Y sin- 
glet and the CEF singlet states, according to a series 
of experiments of linix~o C(T)/T for systems of solid 
solution, Ui-^T^Beia, where the lattice constant ao is 
changed in a wide range covering both the K-Y singlet 
and the CEF singlet ground states. 15 ) Moreover, pure 

UBeis exhibits the NFL behavior, C(T)/T r log T 

up to H z = 12 Tesla. 14 ) Of course, precisely speaking, 
results of the present paper are for the system of f 2 - 
impurity so that we should be careful in deriving a solid 
conclusion. Indeed, an approach based on the dynam- 
ical mean field concept is indispensable for deriving a 
solid conclusion for lattice systems, in which the present 
results would be inherited to the solver of impurity prob- 
lem. Nevertheless, we expect that the effect of the compe- 
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tition plays an important role for UBei3 to exhibit such 
a NFL behavior rather robust against the magnetic field 
larger than the effective Fermi energy inferred from the 
value of liniT~o C(T)/T. Namely, the lower Kondo tem- 
perature would be larger than 12K from the fact that 

the NFL behavior lim T ^ C(T)/T r log T in UBei 3 

is robust against the magnetic field up to 12 Tesla at 
least. 14 ) Were it not for the superconducting state at 
T < T c ~ IK, there would exist the peak with specific 
heat near at T = Tp. Predictions of the present paper 
may be checked by experiments in some U-diluted system 
of UBei3 near the phase boundary between the K-Y sin- 
glet and the CEF singlet states under pressures and/or 
magnetic fields. 
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Appendix: f 2 States in Tetragonal Symmetry 

In the tetragonal symmetry, wave functions for each 
CEF level are given within f 1 states of j = 5/2 orbitals 
as follows: 



.(i) 



) = CI ± |> ^ |> 



irfi) = 



|r 6 ,±> 



±i7|±|)tCIt|), (A-2) 

l±5>, (A-3) 

where £ and r] are the coefficients determined by the 
effect of the CEF. In this appendix, we determine these 
coefficients on the basis of the condition that the energy 
level of low- lying f 2 states with J = 4 manifold can be re- 
produced by Hamiltonian (9). First, we construct states 
with J = 4 manifold from the direct product of states 
with j = 5/2 manifold. 



±3) 



±2 > = ± ^i ± i>i=F5> ± Vni ± i>i ± 5>' 



(A-4) 
(A-5) 
(A-6) 



By using the inversion relation of eqs.(A-l)-(A-3) repre- 
senting | ± 5/2), | ± 3/2), and | ± 1/2) in terms of r^ 1} , 

r^ 2 \ and T 6 , we obtain the f 2 states in the tetragonal 
symmetry as follows: 



1^ = ^(12) 



1 u 



2» 



9jJ L (sc + Vs„) (|r«>|r 6 _> + |r«>|r 6+ >) 
+ (V5( - 3 V ) (|if2)|r 6+ ) - |r£>|r 6 _>)] , (A-7) 
|r 3 > = -J=(|2> + |-2» 

= ^±= [(3< - V5 V ) (|r«>|r 6 _> - |r«>|r 6+ >) 

+ (V5C + 37,) (|r#>|r 8 _> + |if2>|r 6+ >)] , (A-8) 
l^ 2 i) = /?|3>-a|-i> 



«|r«>|rf_>) 



|lf!>=/?|-3>-a|l> 



|r£>|r fl+ >, 



= - [K + ^a V \ |r«>|r 6 _> + ^a|nv>|r 



(A-9) 



->(2)v 
- 7+1 



(A-10) 



(A-l) 



are the same as those expressed in eqs. (5)- (8). 

Here, terms where states with symmetry are occu- 
pied in eqs. (A- 7)- (A- 10) can be negligible because their 
energy levels are assumed to be higher than the other 
states so that the hybridization between and the in- 
states (A- 7)- (A- 10) may be neglected for forming a heavy 
fermion state as discussed in refs. 16 and 17. The coeffi- 
cients, £ and 77, can be determined by the condition that 
the coefficients of the remaining terms in eqs. (A- 7)- (A-8) 
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are equal to those in eqs.(5)-(6). The result is 




(A-ll) 



A relation between a and f3 is also derived by comparing 
eqs.(A-9)-(A-10) with eqs.(7)-(8) as follows: 

VU V 7 V 14 v 7 

It is noted that the coefficient of the first term in (A- 7), 
including becomes larger than that of the second 

term in (A- 7), including T^, if we use the values of 
(A-ll). However, it is allowable to discard the first term 
because the 1^2 state is assumed to play a negligible 



role in forming the heavy fermion state as discussed 
above. The normalization condition for the right part of 
eqs.(A-9)-(A-10) requires x = 1. Nevertheless, by com- 

(2) 

bining the normalization condition for 1^ ; symmetry in 
f 2 states, i.e., \a\ 2 + \/3\ 2 = 1, there are no solutions for 
these coefficients as far as 2vTl/7 < x < 1 (2VTI/7 ~ 
0.947 • • • ). This is because we have discarded the states 
relating to as discussed above, and increased the 
weight of the remaining terms in eqs.(A-9)-(A-10). In 
view of such a situation, for simplicity, we use the pseu- 
dospin representations, | t>t) an d | written in 

eqs.(7)-(8), respectively, as Ig states instead of using 
a and f3. 



